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Fisher zeros are the zeros of the partition function in the complex )3 = INc/g^ plane. When they 
pinch the real axis, finite size scaling allows one to distinguish between first and second order 
transition and to estimate exponents. On the other hand, a gap signals confinement and the method 
can be used to explore the boundary of the conformal window. We present recent numerical 
results for 2D 0{N) sigma models, 4D U{\) and SU (2) pure gauge and SU (3) gauge theory with 
Nf ~ 4 and 12 flavors. We discuss attempts to understand some of these results using analytical 
methods. We discuss the 2-lattice matching and qualitative aspects of the renormalization group 
(RG) flows in the Migdal-Kadanoff approximation, in particular how RG flows starting at large 
j3 seem to move around regions where bulk transitions occur. We consider the effects of the 
boundary conditions on the nonperturbative part of the average energy and on the Fisher zeros for 
the ID 0(2) model. 
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1. Content of the talk 

With the ongoing effort at the LHC, there has been a renewed interest in the phase diagram of 
lattice gauge theory models. The location of the conformal windows for several families of models 
have triggered intense discussions. Different numerical and analytical techniques have been applied 
to QCD-like models with a large number of fermion flavors |jl|, ^, 0, ^ ^, 0] or with fermions 
in higher representations [||, ^, 10, 11 1. See also [|l^, 13, [l4|, [15]] for recent reviews of results and 



expectations. It is important to understand the critical behavior of lattice models from various points 
of view. It was proposed to consider complex extensions [ |I6| , [Ft] , [T^] of the picture of confinement 
proposed by Tomboulis [|T9|]. It was observed that the Fisher's zeros, the zeros of the partition 
function in the complex jS plane, determine the global properties of the complex RG flows. In the 
case where a phase transition is present, the scaling properties of the zeros |2^, 21, ^ allow us 
to distinguish between a first and second order phase transition. In the following, we briefly review 
the Finite Size Scaling (FSS) of the Fisher's zeros. We then discuss numerical results for U{\) 
in 3 and 4D, SU{1) in 4D with ^Adjoint, SU{3) in 4DNf = 4 and 12 flavors. We also discuss the 
effect of boundary conditions in the 0(2) sigma model and recent effort to relate these questions 
to perturbative expansions. 

2. Fisher's zeros and Finite Size Scaling (FSS) 

Fisher's zeros provide information about FSS. The basic principle is the decomposition of the 



partition function [ 24 1 at finite volume into a singular and a regular part: 



Z = Z,,„ge^'""'"*'' (2.1) 

Zsing. = e-^°^""«' (2.2) 

Under a RG transformation, the lattice spacing a increases by a scale factor b, and 

L L/b (2.3) 

fsing. ^ ^ fsing. (2.4) 

^sing. ^ ^sing. ■ (2.5) 



This has the important consequence [ ]25| ] that the zeros of the partition functions are RG invariant. 
The nonlinear scaling variables (e. g. u = P — j3c, . . .) transform homogeneously: Ui — > For 
the relevant variables, we use the notation A, = Z;^/^', and for the iiTclevant variables Ay = b^^J. 
The RG invariance of Zs,„g means that 

Zsing. = QiluiL'^"'}, {ujL-'"^}) (2.6) 

For a single relevant variable m ~ j3 — jSf, we have Zv,„g = Q{uOl^). The complex equation 
Z = can be written as two real equations for two real variables and generic solutions are isolated 
points. 

Z = 0^mL^/'' =w, with r= 1,2,... (2.7) 
This implies the approximate form for the zeros: 

i8,(L)~j3, + w,L-i/^ (2.8) 
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Figure 1: Fisher's zeros for t/(l) on lattices (L=4 (left) and 6 (right)). The zeros of the real (imaginary) 
part are represented by the blue (red) curves and the region of confidence is below the green line (zeros near 
or above this line are not reliable). 



There are many examples, where these discrete solutions follow approximate lines or lay inside 
cusps. In the infinite volume limit, the set of zeros may (or may not) separate the complex plane 
into two or more regions. For a first order transition V is replaced by 1/D. For confining models 
where there is no phase transition on the real coupling axis between weak and strong coupling, it is 
sometimes possible to construct numerically RG flows in the complex j3 plane [ |T6| , [l7| , |T8|]. These 
examples show that the Fisher's zeros can act as "gates" for the complex RG flows and govern their 
global behavior. Since it is much easier to calculate the Fisher's zeros than the complex RG flows 
we shall now describe model calculations of the zeros. 



3. Model calculations 

We showed previously that for 2D 0{N) models in the large-A'^ limit, the RG flows go directly 



from weak coupling to strong coupling [16, [ITp. This means that the theory has a mass gap. We 
would expect the same result for 3D ?7(1). Consistently, using canonical simulations and histogram 
reweighting, we found no zeros near the real axis. This is illustrated in Fig. |I] where anything above 
the green line is considered nonreliable. 

The case of the 4D U{\) is more delicate. There is a transition but is it first or second order? 
Despite some earlier suggestions of a second order transition [^, 22], there seems to be a general 



agreement that the transition is first order [2^, 27, 28]. This question has been revisited recently 
from the point of view of Fisher zeros [|^]. For lattices, the average plaquette distribution has a 
double peak distribution with equal heights at a pseudo-critical jSj. A double peak does not always 
mean a first order transition. We consider a simple example where the plaquette distribution at jS^ 
is a superposition of two Gaussians: 

«(S)e-fe^oc (e-(i/2<T2)(s-s,)^+e-(i/2^')(^-^^)') (3.1) 

The zeros are located at /3,- = j85 + /27r(2r+ l)/(52 — S\). If {S2 — S\) oc lP^ we have a first order 
phase transition (latent heat) and Im/3i oc L^^ . However, if {S2 — Si)°<. L^^^ , then the width of the 
double peak in the average plaquette goes to zero at infinite volume and ImjSi L^'/^ with v = 
1/(D — For the actual U{X) ^"^^^ foi" small L, the distance between the peaks slowly decreases 
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with increasing volume. In the infinite volume limit, the width of the double peak distribution 
of the average plaquette goes to a nonzero limit (latent heat) for a first order phase transition and 
to zero as an inverse power of L for a second order transition. At this point, better statistics at 
large volumes are necessary in our numerical construction of the density of states to discriminate 
between the two scenarios, however fits of the imaginary part of the lowest zero described in Ref. 
[29] show a goodness of fit Q=0.43 for L^^ leading behavior and Q < 10^^ for L^^"^ which favors 
the first order scenario. 

For pure gauge SU (2), it is possible to create double peak distributions by adding an adjoint 
term with a coefficient ^Adjoint sufficiently large and positive to the Wilson action. As a positive 
pAdjoint is increased, the Fisher's zeros go down as illustrated in Fig. ^. There is a clear change of 
behavior for 1.0 < pAd joint < 1.1. Subleading effects are important and larger volumes are being 
studied in order to draw conclusions regarding the scaling of the real and imaginary part of the 
lowest zeros. 

We have started investigating the possibility of doing 2-lattice matching for SU (2) with a 
mixed action [31]. A first step consists in considering the effect of pAd joint on RG flows starting at 
weak coupling. This question can be addressed using the Migdal-Kadanoff (MK) approximation. 
Using expansions with 20 characters, we found that the RG flows seem to go around possible 
boundaries in order to reach the strong coupling fixed point (all /3's = 0). Projections in two planes 
are shown in Fig. |3| and are in qualitative agreement with this idea and flows described by ||30[]. 

It is also possible to test the accuracy of the MK approximation by checking the 2-lattice 
matching. If we had an exact RG transformation between a fine (2L)'* lattice and a coarse 
lattice, then the 2R x 2R Wilson loops for the (2L)'^ lattice and the R x R Wilson loop on a 
lattice with effective couplings obtained by the RG transformation should match exactly. However, 
if we use the MK approximation, the matching is not very accurate especially at large coupling 
as shown in the Table |l|. This could also be used as a test for the Cheng-Tomboulis improvement 
proposed at this conference [ |32| ] . 

The Fisher's zeros for SU (3) with A'^y= 4 and 12 quarks in the fundamental representation have 
been investigated [33]. We used the standard Wilson action and unimproved staggered fermion 
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Figure 3: Projection of the MK flows in the {^1/2, Pi) and (/33/2:j3i) planes. 
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Table 1: Wilson loops on fine and coarse lattices. 



action with the Rational Hybiid Monte Carlo (RHMC) algorithm. We started with relatively small 
symmetric lattices and up to 50,000 configurations. The bare quark mass is set to be = 0.02 for 
now. For Nf = 12, we found a discontinuity for the plaquette near /3 ~ 4. 1 for V = 8^ lattices as 
shown in Fig. ^ left panel. The plaquette histories showed the characteristic hysteresis behavior. 
For comparison, we show the crossover for A'^^ = 4. The lowest zeros are shown in Fig. ^ right 
panel. Possible hypothesis to be tested are that for Nf = 12, the imaginary part scales like L^^, 
which signals a first order phase transition and that the real part increases like log(L). Nonlinear 
effects seem important at small volume and larger volume calculations are in progress. As already 
noticed in Ref. [|l|], the chiral condensate has a discontinuity near the same value of /3 as the 
plaquette. More recently, it was shown that improved actions can create a second discontinuity 



with a broken single-site shift symmetry between the two transitions [34, 35]. 



We also studied the Fisher's zeros for the ID 0(2) with L = 4, 8, 16, 32 [|36|]. As shown in Fig. 
gthe zeros are very different for open (o.b.c) and periodic boundary conditions (p.b.c). The MK 
complex flows can also be constructed and are also shown in Fig. ^. At finite volume, the nonper- 
turbative parts of the average energy are very different for open and periodic boundary conditions. 
It was found that \{E - Ept)/E \ oc q-^-P for open boundary condition and e for periodic 
boundary condition, where E^. is the energy of the periodic solution of the classical equation of 
motion with winding number 1. Hadamard series have been constructed to improve the accuracy 
at strong coupling. 
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Figure 4: Average plaquette and Fisher's zeros for SU (3) with Nf= 4 and 12 at different volumes. 




Figure 5: MK complex flows for the system (o.b.c) and zeros from the two different boundary conditions 

4. Conclusions 

Much progress has been made in finding reliable ways to locate the Fisher's zeros of various 
models. A consistent picture of confinement in terms of complex RG flow is emerging. Much work 
remains to be done for understanding multiflavor models within this approach. Better analytical 
approaches (based on improved RG or weak coupling expansions) are needed. The FSS of zeros 
is simple, however subleading connections are important (at least for unimproved actions). We plan 
to monitor the effects of improvement on the zeros by turning on improvement adiabatically. 

This research was supported in part by the Department of Energy under Contracts No. 
FG02-91ER40664, DE-FG02-97ER41022, DE-FG02-12ER41871 and DE-AC02-98CH10886. 
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